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Abstract 

The recently proposed states of low energy provide a well-motivated 
class of reference states for the quantized linear scalar field on cosmological 
Friedmann-Robertson-Walker spacetimes. The low energy property of a 
state is localized close to some value of the cosmological time coordinate. 
We present calculations of the relative cosmological particle production 
between a state of low energy at early time and another such state at 
later time. In an exponentially expanding Universe, we find that the par- 
ticle production shows oscillations in the spatial frequency modes. The 
QQ ' basis of the method for calculating the relative particle production is com- 

pletely rigorous. Approximations are only used at the level of numerical 
calculation. 
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tl' 1 Introduction 

The phenomenon of particle production on non-stationary cosmological space- 
times is one of the early findings in the theory of quantized fields on curved 
^ , spacetimes [T7]. At its very basics, it can be reduced to an analogy between 

the equation of motion for the mode functions of a quantum field in a time- 
, dependent (spatially homogeneous) spacetime metric and the equation of motion 

for a harmonic oscillator with time-dependent frequency [17118]. It is, however, 
, not easy to make the analogy completely rigorous, for the mode decomposition 

' on which the notion of particles is based is itself time-dependent. 

0^ , Quite generally, in a time-dependent background, any concept of particle 

' will itself depend on time, and possibly on other data, such as scale parameters. 

In more precise terms, for a quantum field on a spacetime with time dependent 
metric, a concept of particle is tied to the selection of a reference state, which 
then depends on such data as a hypersurface of constant time, and other, soak- 
ed ' setting parameters. In the present paper, we will be concerned with the case of a 
linear, minimally coupled scalar field on Friedmann-Robertson-Walker (FRW) 
spacetimes. For this situation, particle creation with respect to some classes 
of reference states has been investigated. Notably, Parker introduced so-called 
"adiabatic vacuum states" , designed to minimize particle creation, as reference 
states [T71 H]. The precise mathematical definition of this class of states is 
intricate and was first achieved in an investigation by Liiders and Roberts [15j . 
In further work by Junker and Schrohe JjJj , it was demonstrated that adiabatic 
vacua are locally unitarily equivalent to the class of Hadamard states. Since 
many results have, over the past decades, confirmed the view originally put 
forward by Wald [531 US] that the class of physical states for linear quantum 
fields on curved spacetimes should contain all (quasifree) Hadamard states as 
building blocks, this shows that adiabatic vacua qualify as physical states as 
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well. However, the physical interpretation of adiabatic vacua remains less clear; 
furthermore, their definition involves, in the general case, asymptotic series 
expansions making them inconvenient to handle and posing difficulties when 
trying to obtain numerical results on particle creation |27j . 

Quite recently, Olbermann [1^ has introduced a new class of reference states 
for the linear scalar field on FRW spacetimes, called states of low energy (SLEs). 
These states are spatially homogeneous (with respect to the spatial isometry 
group of FRW spacetimes) and such that that they minimize the time-integral 
of the energy-density for a given weighting function. For the moment, let us 
describe this in more detail as follows (the exact definition will be given in Sec. 
3). Suppose that t is the cosmological time coordinate in an FRW spacetime, 
and let f{t) denote a non- negative, compactly supported, smooth function of 
cosmological time. For a homogeneous state of the quantized linear scalar field 
on FRW spacetime, denoted by its expectation functional ( . where w is a label 
for the state, we write (: Q(t) for the (renormalized) expected energy density 
in that state. As the state is homogeneous, this quantity is only dependent on 
time t. The time-integral of the expected renormalized energy density, weighted 
with respect to /, is 



Then a state ( . is called an SLE (state of low energy) with respect to the 
weighting function f if the state minimizes the expression g^^ [f] among all ho- 
mogeneous states {■)uj- Note that the SLE property depends on the choice of 
weighting function / and that it is therefore a concept that is local in time. Ol- 
bermann has shown that SLEs are Hadamard states, and that adiabatic vacua 
are approximations of SLEs. 

If a state {.)^ were an SLE for the case / = 6tg, the delta-distribution 
concentrated at cosmological time t, it would correspond to a state minimizing 
the (homogeneous) energy density at time to. This is reminiscent of the concept 
of an instantaneous vacuum state at time to, and the interpretation of SLEs is 
that of approximate (or, by weighting with respect to /, mollified) instantaneous 
vacua at time to- Instantaneous vacua, however, aren't Hadamard states, and 
they are not even (locally) unitarily equivalent to Hadamard states [Q] [12], 
thus disqualifying them from the class of physical states. The weighting, or 
mollifying, of the energy density by integrating in time against smooth weighting 
functions / in defining SLEs is therefore needed in order to ensure that SLEs 
are physical states. 

Nevertheless, the interpretation as homogeneous, energy minimizing states 
close to an instant of cosmological time makes SLEs natural candidates for ref- 
erence states with respect to which particle concepts can be defined and particle 
creation in FRW spacetimes can be calculated. This is the topic of the present 
work. More specifically, we consider two smooth weighting functions, /i and 
/2, peaked at cosmological time parameters ti and t2, respectively. We will 
envisage the situation that ti corresponds to an early cosmological time and 
t2 > ti a later cosmological time. Denoting by and {■)uj{2) the cor- 

responding SLEs, we will determine the particle content per (fixed) frequency 
mode k (arising from the spatial symmetry) of the states with respect to each 
other; this quantity is given by the modulus squared |/3fcp of the Bogoliubov 
coefficients, see Sec. 4. The numerical calculation, however, turns out to be 
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involved, and therefore we limit ourselves in this work to a special case which 
simplifies the calculations significantly. However, we think that the principal 
qualitative aspects of the results of our calculations won't depend very critically 
on these choices. Our specialized assumptions are: (i) A closed (spatially com- 
pact) FRW spacetime, (ii) an exponentially increasing scale-factor a{t) — e^* 
{H > 0), (iii) a very small mass parameter m > 0, fine-tuned with respect to 
the Hubble-parameter H, in the linear field equation of the quantized scalar 
field. Moreover, we approximate the weighting functions /i and /2 by Gaus- 
sians peaked at ti and t2 with small characteristic widths e and 5, respectively. 
Then will be calculated numerically and represented graphically for certain 
choices of ti,t2,e and S. 

The most interesting result is that |/3fcp shows oscillations in k when t2 — ti 
is different from and this effect increases with growing ^2 A somewhat re- 
lated observation to the effect that the particle number per mode can decrease in 
time appears in llj however, the reference states used in [llj are instantaneous 
vacua. 

It is worth pointing out that our appraoch is completely rigorous, with ap- 
proximations entering only at the level of numerical calculation. While the 
simplifications (i), (ii) and (iii) made in order to facilitate the numerical calcu- 
lations are not physically realistic, it is, as mentioned, to be expected that the 
basic findings of our calculations, especially the said oscillatory behaviour of 
will qualitatively also occur in physically more realistic situations. This 
could have some implications with regard to observations in cosmology. 

This article is organized as follows. The quantization of the linear (minimally 
coupled) scalar field on FRW spacetimes will be summarized in Sec. 2, together 
with the definition of the renormalized stress-energy tensor. States of low energy 
will be discussed in Sec. 3. In Sec. 4 we consider Bogoliubov transformations 
between two SLEs and the associated notion of particle creation. The numerical 
calculations will be presented in Sec. 5. Summary and Outlook are given in the 
final Sec. 6. 

2 Quantum field on FRW spacetimes and renor- 
malized stress-energy tensor 

In this section, we summarize the quantization of the linear scalar Klein-Gordon 
field on FRW spacetimes, and the definition of the renormalized stress-energy 
tensor. This serves mainly to make the present text as self-contained as possible. 
General references for the material in this section are [531 IMl US [S] . 

A spacetime {M,g), where M denotes the spacetime manifold and g the 
Lorentzian metric, is of FRW type if it is of the form 

1. A/ = R X E 

2. dsl = dt^ - a{tfK.jdy.'dyJ 

3. (S, h) is a 3-dimensional manifold with Riemannian metric of constant 
curvature which is either equal to or normalized to ±f. (We write 
usually X for elements in M, identified also as (i,x) with t e M, x e S.) 

4. a : M ^ R is a C°° function taking strictly positive values, called the scale 
factor. 
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The classical linear, minimally coupled scalar field (Klein-Gordon field) on 
{M, g) fulfills the field equation 

(□g +m^)^ = 

with smooth ip : M 'R. The constant parameter m > is called the mass 
term of the field equation and Dg is the d'Alembertian of (A/, g) which in the 
present case takes the form 

□.+n.^^g + 3F(t)|--i^A. + .n^ (1) 

Here, H{t) = {-^a{t)) / a{t) is the Hubble function, and Ah denotes the Laplace 
operator on (S, h). An FRW spacetime {M, g) is a globally hyperbolic spacetime 
and therefore there exists a unique pair, E^, of advanced/retarded Green's 
functions for the Klein-Gordon operator Dg -I- m? , for any fixed m > [SI [I]- 
They are characterized as continuous linear maps : (M, M) — > C°° (Af, M) 
with E^{Ug + to2)f = f = (Dg m^)E^F for aU F e C^(Af,R), and with 
supp(£;=^F) C J=^(suppF) for all F e C^(Af,R), where J^(G) denotes the 
future(-|-)/past(— ) causal set of a subset G C M. The difference E = E^ — E^ 
is often called causal propagator of the linear scalar field on (Af, g) . We write 

£(Fi,F2)- / Fl{x){EF2){x)d^ig{x) (Fi,F2eCo°°(M,M)) 
JM 



for the associated bilinear form, where d^g{x) — \dei{g fiu)\d'^ x is the metric- 
induced volume form on M . One can show that £ is antisymmetric, i.e. 

£(Fi,F2) = -£(F2,Fi). 

The linear scalar field on {M,g) can then be quantized as follows. One 
defines an abstract *-algebra ?(Af, g) with algebraic unit element 1 as being 
generated by 1 and a family of elements 0(f), f G {M,M.), with the relations 

0(Fi)(/>(F2) - <?!)(F2)0(Fi) = f£(Fi,F2), 

(j){{ag +m^)F) = 
for aU F,Fi,F2 e Co°°(Af,R). 

One can show that such an algebra exists (nontrivially) and that it is unique 
up to (natural) isomorphisms (cf. e.g. and lit. cited there). The generating 
elements 0(f) can be viewed as "abstract field operators" for the quantized 
linear scalar field on {M,g). The (j){F) can be turned into operators acting in a 
Hilbert space upon considering Hilbert space representations of 3^{M, g) . Here, 
we are interested in Hilbert space representations induced by states on 3^{M,g). 
Recall that a state on 3^{M, g) is a linear functional w : 3'(Af , g) C which 
is positve, i.e. it fulfills uj{X*X) > for all X € 'J{M,g). (Note that X is in 
?(M, g) if it is a polynomial built out of 1 and finitely many 0(fi), . . . , (^{fn)-) 
Moreover, it is required that u) be continuous, that is, for each n e N the map 
Fi (g) • • • (g) F„ 1-^ u;(0(fi) • • • 0(f„)) extends to a distribution in (C(f'(Af R))'. 
It is also common to denote a state by its expectation value functional, i^X^^^ — 
[x'(X), and we shall often adopt this notation. 

There are states u entirely determined by their 2-point function 

W^(Fi,F2) =L^(0(Fi)0(F2)) = (0(Fi)0(F2))^ (Fi,F2 G {M ,^)) 
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by requiring the relations (0(fi) • • • 4>{F2n+i))Lu = and 

for all n e N; the left hand side of this equation is to be read as the expectation 
value of the polynomial in (/i(f) that results from formally differentiating e*'^'*'^^' 
and setting A equal to 0. 

A state u! on 3^{M, g) induces a ^-representation tt^^ of 3^{M, g) on a dense do- 
main in a Hilbert space "K^ together with a canonical unit vector € . 
The collection of {'Kun'O^.TT^,^^^) is called GNS-representation, or Wightman- 
representation of uj (cf. [H]) and is characterized by the properties that = 
{n^,n^{X)^l^) for all X G J(M,<7), tt^{X*)x = 7r<.(^)*X for all x £ 2^^, and 

— {i^ujiX)^^ : X G 3'(M, 5)}. Then the represented abstract field operators 
(/)(/) become unbounded operators 

^M) = -^^im) (2) 

defined on the domain in the representation Hilbert space yiui- 

For a quasifree state uj, there is a one-particle Hilbert space JCif'' together 
with a real-linear map : Cg°{M,R) Jfi^^ so that {K^{fi), K^{f2)) = 
W2(fi,F2). One can then show that Jf^j = F+(!K£^''), the bosonic Fock-space 
over the one-particle space JCii^^ , and that 

7r^(0(F)) = a(if^(F)) + a+(K^(F)) (3) 

where a{.) and a*(.) are the usual annihilation and creation operators on 
Furthermore, = (1,0,0, . . .) is the Fock -vacuum- vector and thus 
it is convenient to write \0)u, — for the GNS-vector of a quasifree states uj. 
Moreover, a quasifree state uj is pure if the range of is dense in This is 

equivalent to saying that uj cannot be written as convex combination of several 
different states. 

On each FRW spacetime {M = M x S,^) there acts the spatial symmetry 
group Gs- It acts only in the S-part of M, so that for each 7 G Gs one has 
7(t,x) = (t,7x;(x)) for all i G M, x G E with respect to the R x S-splitting 
of M, with an isometry 75] of (E, By setting a^((/)(F)) = 0(f o 7^^) for 
F G C^{M,M.) and 7 G Gs, one can define the automorphisms : S'{M,g) 
3^{M,g) induced by 7 on the algebra of abstract field operators. Then a state uj 
on 3^[M,g) is called homogeneous if it is invariant under the action of a-y, i.e. if 

{a^{X))^ = (X)^ (7 G Gs, X G HM,g)) . 

If w is a quasifree state, then it is homogeneous exactly if 

W^(Fi o 7-1, F2 o 7-I) = W^(Fi, F2) 

holds for all 7 G Gs and Fi, F2 G G^(M,R). 

The most important class of states for quantized linear fields on curved 
spacetimes, particularly in the context of defining expectation values of the 

^Strictly speaking, Ge = {idiR X 7E • 7E G Iso^(E,/i)} where Iso^(S,/i) is the Lie-group 
of orientation-preserving isometries of (S,/i). 

5 



,-AW-(F,F)/2 (f e C'tf'CM,]!)) 



A=n 



stress-energy tensor, are quasifree Hadamard states. These are quasifree states 
whose 2-point functions are of so-called Hadamard form. We shall not go into 
full details of the definition of Hadamard form (see [T11[5S1[7]) for further dis- 
cussion) but mainly describe its basic features entering into the definition of the 
renormalized expected stress-energy tensor. Basically, is of Hadamard form 
if 

W^(Fi,F2) = gsing(Fl,F2) + J J Fi{x)R^{x, x')F2ix') d^Ig{x) dfig^x') 

where Ssing is a singular contribution (a distribution in (C^(M x M))') which 
depends in a certain, local way on the spacetimc metric g and the mass param- 
eter m in the Klein-Gordon equation, but is independent of uj (i.e. the singular 
part is the same for all Hadamard states lo). There remains a smooth contri- 
bution, expressed by Ru, G C°°{M x M,C), which contains the dependence on 
the states u. The circumstance that the singular part 9sing is the same for all 
Hadamard states lo is instrumental for the definition of the expectation value 
(: Tpjy :)(^ of the renormalized stress-energy tensor in the state uj. 

We shall very briefly elaborate on the "symmetric Hadamard parametrix" 
(SHP) renormalization of the expectation value of the stress-energy tensor which 
was employed in |19| , see also [7] . Define the symmetric Hadamard parametrix 

1 i 

S(Fi,F2) = -(gsing(Fl,F2) + Ssing (F2 , F2 ) ) + -£(Fi,F2) 

and, correspondingly, set 

3iE"''(Fl,F2)=W-(Fi,F2)-g(Fi,F2) 

for Fi,F2 e Co°°(M,M). Then the distribution Ji^^^ is given by a smooth 
integral kernel 

3?5HP(Fi,F2)= J jF,{y)Rl^^{x,x')F2{x')dflg{x)df^g{x'). 

Now let X S AI be given and suppose that x' G M lies in a convex normal neigh- 
bourhood of X, and denote by y v^x'; x) the operation of parallelly transporting 
a covector f^, at x' to a covector at x. Note that ^{x'; x) depends smoothly 
on X and x' . With this convention, one can define the SHP-renormalized expec- 
tation value of the stress-energy tensor in the state oj as 

(TSHP(x))^ = hm (V^V.,i?5HP(x,x'))r'''.(:r';x) 

+ ^.g^''(a;)(V^V.,i?EHP(x, x'W'Ax'; x) 
+ ii?5"P(x,.'). 

Note that on the right hand side, operates with respect to x and operates 
with respect to x' . The resulting (T^^^(a;));^ is a smooth, symmetric co-tensor 
field of cc € M . This renormalization procedure has the advantage of being 
completely independent of any "reference" state; in fact, it renders a local, 
generally covariant quantity in the sense of [251[5J[Tn]. However, in general it 
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will have the defect of not being divergence- free. As pointed out in [231 US], this 
defect can be repaired as follows: One can show that there is a smooth function 
Q on M, determined entirely by the local geometry of {M,g), such that 

V''(TfHP(a:))^ = V.Q(:r) {xeM) 

holds for all quasifree Hadamard states oj - while Q is state- independent. Hence, 
if Q{x)gf_ii,{x) is subtracted from {T^^{x))i^, the resulting quantity has vanish- 
ing divergence. In fact, this resulting quantity is a local, generally covariant and 
divergence-free definition of the expectation value of the stress-energy tensor of 
a quantized field which thus complies with the requirements delineated by Wald 
for stress-energy observables of quantum fields [231 ES] ■ Then there remains a 
renormalization ambiguity for this quantity to the effect that one may add other 
symmetric co-tensor fields C^j^ having vanishing divergence and determined by 
the local geometry of (M, g) . Following [7] , we adopt the point of view that the 
specification of C^i^ is a further datum (akin to the mass parameter m) of the 
quantum field (p. Using an argument of [25], the freedom of choosing C^i^ can 
be further narrowed down, so that 

C^. = C^AA r, D] = Ag^, + BG^, + ^j^Si{g) + D^S^ig) , 

where G^i, is the Einstein tensor, Si{g) = j^jR^d^ig, S2{g) = Jj^j R^i.R^^" dfig 
with R and R^^, denoting scalar curvature and Ricci tensor, while 5/5g^'^ denotes 
functional differentiation; the real constants A, _B,r and D parametrize the 
remaining renormalization ambiguity. In conclusion, the renormalized expected 
stress-energy tensor in a quasifree Hadamard state to takes the form 

(: T^, : (x))^ = (T^HP(x))^ - Q{x)g^,{x) + C^^M. B,T , D]{x) {x G M) . 

Considering the case of an FRW spacetime with cosmological time coordinate 
t, the renormalized expected energy density in a quasifree Hadamard state ld 
with respect to cosmological time is 

(: Q : (t, x))^ = (: T : (t, x))^ (^) (|) ((^^ x) G R x S) . 

It is worth making a few remarks here. 

(i) The notation using double dots is to be understood as signifying that 
(: T fj,u '■ is a renormalized quantity. It is not to be confused with the more 
common usage of indicating "normal ordering" , which refers to a reference state. 
The same applies to our notation of (: g : {x))i^ for the renormalized expected 
energy density. 

(ii) The above indicated procedure of defining the renormalized stress-energy 
expectation value of a linear quantized field in Hadamard states u applies not 
only to the case of FRW spacetimes but to general, globally hyperbolic space- 
times. 

(iii) It follows from the appearance of renormalization ambiguities in the def- 
inition of (: T^i, : (parametrized by renormalization constants A, B, F, D) 
that there is no intrinsic (without further considerations) prediction of the total 
absolute value of the (expected) local energy density in quantum field theory 
on generic spacetimes prior to fixing the renormalization ambiguities. The way 
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in which the fixing is done can have significant impHcations in cosmological sce- 
narios (see, e.g., [3J). In the context of discussions of the role played by the 
"vacuum energy" as contribution to the cosmological constant it is occasionally 
claimed that quantum field theory were predicting a value of the vacuum en- 
ergy which misses the observed value by 120 orders of magnitude. It is worth 
pointing out that quantum field theory doesn't make any such prediction. 

3 States of low energy (SLEs) 

Since our investigation later in this article refers to the case of spatially compact 
FRW spacetime, we will from now on restrict our discussion to that case, mostly 
to simplify notation. However, states of low energy have been introduced in [16] 
for all types of FRW spacetimes. 

The spatially compact FRW spacetime has Y, = with the Riemannian 
metric on E derived from the embedding 

i=i 

into equipped with the Euclidean metric. This Riemannian manifold carries 
the metric-induced measure dfih which in turn gives rise to the Hilbert space 
L^{T,,diih)- The Laplacian Ah is essentially selfadjoint on C°°(E, C) and there 
is an orthonormal basis {Ik}, k = {k,l,m), k G Nq, I — 0,1,..., fc, m = 
— /, — Z + 1, . . . , Z, of C°° functions on E which are eigenvectors of — A/j, 

- A,,Yk = At(k)2rk , K{k) = ^k{k + 2) . 

On spatially compact FRW spacetimes, all homogeneous pure quasifree states 
can be represented in a particular form, and this will be the starting point for 
the definition of SLEs. To this end, we quote the following result. 

3.1 Theorem. (Liiders and Roberts |15j . Olbermann [16]) 

(A) Suppose that uj is a homogeneous pure quasifree state on ?(Af, g) for a spa- 
tially compact FRW spacetime with scale factor a{t) and Hubble function 
H{t). 

Then there exists a sequence {Tfc}, k € Nq, of C°° functions : R — > C 
which are solutions of the differential equations (where a dot indicates 
differentiation with respect to t and overlining means complex conjugation ) 

fk{t) + 3Hit)fk{t)+Lol{t)Tk{t)^0 (teR) (4) 
with the time- dependent frequencies 

Lul{t)^K{k)/a'^{t)+m'^ (5) 
and which fulfill the condition 

Tk {t)fk it)~tk {t)Tk (t) ^ la"" (t) (t e K) , (6) 
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such that, in the sense of distributions, one has 

w-((t,x),(t',x'))= Yl rfe(t)yk(x)rfc(OFk(x'). (7) 

k— (/c,/,m) 

Relation ([7|) can equivalently be expressed by stating that the one-particle 
real-linear map given as follows: 

^^oj^ — i'^ = space of square- summable sequences {sk} indexed by the 
k — (fc, I, to), 

K^{f) = {K^{F)i,} with 

K^f)^^ J^J^Tk{t)Y,,i^)Fit,^)dtdM^). (8) 

(B) Conversely, let {Tk}, fc G No fee a sequence of smooth functions : M — )■ C 
fulfilling ^ and ^ . Then the right hand side of ([7]) defines the two-point 
function W2 of a homogeneous pure quasifree state to on 3^{M, g) for a 
spatially compact FRW spacetime with scale factor a{t). 

Thus, there is a one-to-one correspondence between the set of homogeneous 
pure quasifree states cu and the set of sequences {T^} of smooth functions ful- 
fiUing dH) and ©H 

Let us, from now on, denote this set of sequences by We shall next 
quote a result by Olbermann which ultimately introduces states of low energy. 
These make reference to weighting functions. We shall say that that a function 
/ : R ^ M is a weighting function if there is a function q € C^(M,R) such that 
fit) = q{tf {t e R). 

3.2 Theorem. (Olbermann [16] ) 

For any weighting function f and any {Tk} € define the sequence of numbers 
r,{f,{Tk}) = Jjitmit)]' +u;^{t)m)\^)dt ij e No). (9) 

o 

(i) Suppose that, for fixed weighting function f , the sequence {Tk} G -3^ has 
the minimizing property 

r,{f, {Tk}) = min r,(/, {T,}) (10) 

for all J G No- Then the corresponding homogeneous pure quasifree state 
U! — LU f on 3^{M, g) is a Hadamard state. Furthermore, this homogeneous 
pure quasifree Hadamard state minimizes the energy density weighted by 
f along cosmological time, i.e. 

ffoj/] = min (11) 

where the minimum is taken over all homogeneous pure quasifree Hadamard 
states Lo, with 

qM]^ I fit){:Q:it,^))u.dt (12) 

^States which here are called "homogeneous" are called "isotropic" in 1161 . 
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for any x G S. (Owing to homogenity of u, the quantity on the right hand 
side of (|12p is independent of x.) 

o 

(a) For given weighting function f , a minimizing sequence {T^} G £^ exits. 
This sequence is unique up to choice of a phase factor for each k, i.e. 

o o 

{Tfc} has the minimizing property exactly if this is the case for {e'^'^Tfe} 
for any sequence of real numbers {pk}- All choices of phases lead to the 

o 

same minimizing state uj. 

o 

Given any sequence {Sk} G a minimizing sequence {Tk} can be con- 
structed from {Sk} in the following way: One sets 

Tk{t) = \(k)Sk{t) + fi{k)Sk{t) , (13) 

with the definitions 



m 



ci(fc) 1 



2Vci(fc)2-|c2(fc)|2 2 
V2Vci(fc)2-|c2(fc)|2 2 



ci{k) = l I f{t){\Sk{t)\'+iul{t)\Sk{t)\')dt 

C2{k) = \j f{t){Sk{tf+ui{t)Sk{tr)dt 

Cfe = TT - ArgC2(fc) 



(14) 



Remarks 

(a) Note that {T^}, respectively to, minimize Qui\f] no matter how the renor- 
mahzation constants A,B,T,D have been chosen. 

(/3) For fixed x G S, /jj /(t) (: q : {t, x.))^ dt equals the /-weighted integral of the 
expected energy density along the geodesic t i— > (t, x) in FRW spacetime. It is 
known that this weighted integral is bounded below as a functional on quasifree 
Hadamard states lo. This fact is a special case of a quantum energy inequality, 
which has been established for the minimally coupled scalar field by Fewster [H] . 

o 

(7) Note also that A(fc) and fi{k) and hence the minimizing sequence {Tk} are 
unchanged under a constant rescaling f(t) 1— > r • f(t) (r > 0) of the amplitude 
of the weighting function. 

4 Bogoliubov transformations and particle cre- 
ation 

We continue to consider the case of spatially compact FRW spacetime (M, g) 
with scale factor a(t) and Hubble function H{t). 

For any pure homogeneous quasifree state w on 3^(Af, g), the one-particle 
map : C^(M,R) -> :k£^^ has Jf£^^ = with K^{f) given by ©; is 
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identified as the space of square summable sequences {sk}, k = {k,l,m). The 
field operators <&i^(f) — 7r^(0(F)) in the GNS-representation of uj take the form 

$„(f) = a(if^(F))+a+(/f^(F)) (15) 

by and ([3]). Let us denote by {(5k', k} the sequence in £^ which takes the 
value 1 exactly if k = k', and the value else. Then write a^i = a(<5k',k) and 
a^, = a^((5k',k)- With this notation, and writing more simply k instead of k', 
one may recast (jisp in the form 

$^(t,x)= (rfc(i)Fk(x)a++Tfe(t)yk(x)ak), (16) 

k— (fc,/,m) 

to be interpreted as an operator-valued distribution (the right hand side be- 
comes an operator in Fock space F_|_(lKi^-') upon integrating with a test-function 

F(t,x)). 

It is known that for any pair of pure quasifree Hadamard states w(l) and 
Lj(2) on spatially compact FRW spacetime there is a unitary operator U : 

F+i^lt)) ^ ^+(^^2)) s'ich that 

^^i2){F)^U<P^^2){F)U-' (17) 

for aU F e Cg°{M,R) [HI HI]. 

Now let a;(l) and U!{2) be two pure homogeneous quasifree Hadamard states 
and denote by {T^^''} and {Tj^^^} the corresponding sequences in Since T^^"* 
and T^^-* both fulfill the differential equation ^ and the solution space of (UJ 
is two-dimensional, there are complex coefficients ak and (3k so that 

The requirement ^ implies 

|afc|'-|/3fc|'-l. 
Insertion in to (I16p then yields the relations 

$^(i)(i,x) = J2 (Ti'\t)Fk(x)a++Ti'\t)Fk(x)ak) and 

k— (A; J,m) 

*.(2)(i,x) - (rf (t)Fk(x)6++Tf (t)yk(x)6k) 

k— (A:,/,m) 

with 

^k = afcoJ + Z^kO-k, (18) 

6k = akOk + Pka\ , (19) 

where the notation — k = (fc, I, —m) for k = (k, I, m) and the property = Y-u 
have been used. The passage from annihilation and creation operators ak and 
Q-k and by (llSp , (I19p is a special case of a Bogoliuhov transformation. 

It is related to the unitary U appearing in (|17p as formulated in the following, 
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well-known theorem ^01 [51 [T5] : 

There is a unitary U : F+i^ij^l^^) ( Jf^^^^^j ) such that 

$c.(2)(f) = (7$„(i)(f)[/-i (f e Co°°(M,R)) 

exactly if 

bk = UauU^^ and b+ = Ua+U^^ 
for all k. This, in turn, holds if and only if 

I](fc' + 2fc + l)|/3fe|2<cx). (20) 

W=(k,l,m) feeNo 

Now suppose that we consider the GNS-representation of a;(2) as "reference" 
representation. With respect to this representation, we have 6k|0)t^(2) = 0, so 
the 6k and 6j are annihilation and creation operators of the spatial symmetry 
mode k with respect to this representation. One refers to this property also by 
saying that 6k and annihilate and create "b-particles" . Correspondingly, the 
number operator 

counts the number of of k-modes with respect to having chosen the GNS- 
representation of w(2) as reference representation. One abbreviates that by 

saying that ^k^^' counts the number of k-modes of b-particles. Thus, J^k^^^ 

be regarded as a counter for k-modes which has been calibrated to give zero 

response in the reference state {■)u>{2)- 

Then rik*'^'' gives |/3/cP as response in the state ( . )ui{2)' 

K^'^)c.(i) = .(i)(0|6+6k|0)^(i) 

= a;(i)(0|(afcak + /3fca-k)(afcak + f3ka\)\0)u,{i) 
= \Pk\' 

observing that au\Q)u,(i) = and a;(i) (0|a^k»-k|0)c^(i) = 1- 

We remark that in a situation where the scale factor a{t) is never constant, 
it is impossible to have |/9fcp = for all k unless the weighting functions /i 
and /2 of the two SLEs w(l) and w(2) are proportional {fi{t) = r ■ f2{t) with 
a positive constant r). To see this, note that /3fe = imphes that 6k|0);^(i) = 
which, in turn, means that ak.U~^\0)i^(i) = 0. But the only non-zero vectors 

in F^{'K''^^-^s^) which are annihilated by all Ok are scalar multiples of |0)^(i). 
Likewise, the only non-zero vectors in (IK^^^2) ) which are annihilated by all 
the 6k are scalar multiples of 10)^^(2). Hence, J7|0)„(i) = e*'"|0)(^(2) for some real 
r. By the properties of the GNS-representation, this implies that uj{l) = uj{2) 
as states on J'{M,g). However, Olbermann (TB] has shown that for this to hold 
with non-constant a{t), it is necessary that the weighting functions /i and /2 are 
proportional. This feature for non-constant a{t) is significantly different from 
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the case of constant a{t) = ao, where there is one unique SLE for all weighting 

functions and hence always /3fc = 0. 

In the case of non-constant a{t), the deviation of Pk from is interpreted 
as (cosmological) particle creation. One may envisage the following situation 
portrayed in Figure 1: Some "initial" state uj{l) is prepared as SLE with respect 
to a weighting function /i concentrated near an early cosmological time ti. 
Another SLE uj{2) corresponding to a weighting function /2 concentrated at 
later cosmological time t2 > ti is used as reference state. This reference state 
can be taken as an approximate vacuum for an observer making measurements 
near the time t2- If the observer uses a particle (mode) counter calibrated to 
give zero response in his "approximate vacuum" i-o{2), he or she will find that 
the symmetry mode k measured on the state u;(l) is excited with a distribution 



On the other hand, another observer making measurements around the time 
ti using a;(l) as reference state will find no particle excitations when making 
measurements on the state The way in which this "particle creation" 

depends on reference states that are defined with respect to properties localized 
in time should be noticed here. 



Figure 1. Illustration of the localizations of the weighting functions fi and 
/2 in an exponentially expanding Universe. The parameters 6 and e denote 
the characteristic widths of these functions. The curve 7 depicts the timelike 
geodesic of an inertial observer. 
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5 Calculation of the number of created particles 



As described in the previous section, we want to calculate the number of created 
particles per mode k. We will make some special assumptions with regard to 
the Hubble function H{t) and the mass parameter m and some further approx- 
imations for the sole purpose of facilitating the numerical calculations. Let us 
explain how we proceed in several steps. 

(I) 

Supposing that a pair of weighting functions /i and /2 has been chosen, we de- 
note by uj{l) and uj(2) the corresponding SLEs (pure, homogeneous, quasifree) 
and by {T^^-*} and {T^^-*} the corresponding sequences in They are related 
by the Bogoliubov coefficients ak and (3k (fc G No) via 

Tf ) (t) = afcTf ^ (t) + PkTi''' {t) {teR,keNo) (21) 
Using moreover (cf. ^) 

T^\t)tl'\t)-TkTi'\t) = ia{t)-' it€R,k€ No) (22) 

where a{t) is the scale factor of the underlying FRW spacetime, one finds that 
(PT|) and (HH) imply 



/5. = -^(T^fr«-r«rf)[, (23) 



for any choice of t* G M. 
(11) 

The best way to find {T^^-*} and {T"^^"*} is to choose some {Sk} G T and to 
calculate the coefficients X^^^k) and fi'-^^{k) {j = 1,2) as in (...) so that 

ri'' it) = A(^') ik)Sk it) + Ai(^) ik)Sk it) . (24) 

Observing that 

SkSkit) ~ '^kit)Skit) = ia{t)-' , (25) 
equations and combine to yield 

Pk = A(i)(A)/x(2)(fc) - A(2)(fc)Ai(i)(A) (26) 

for all fc e Nq. 
(Ill) 

We will now specify the Hubble function as H{t) ~ H with a constant H > 0, 
thus 

ait) = e^' , (27) 

and we will choose the unit of time such that H — 1. In view of the fact that 
we have assumed our underlying FRW spacetime to be spatially closed, this 
scenario doesn't comply with a solution of the Einstein equations with normal 
matter (ideal fluid) without cosmological constant; rather it models an epoch of 
accelerated expansion of the Universe that may have taken place over a time- 
span such that H — 1. For the purposes of illustration, we will nevertheless take 
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that time-span here to be the estimated age of the Universe, tn — 1-39 • lO^^y. 
(One may take any other (shorter) time-span upon scahng the mass parameter m 
accordingly to interpret our numerical results analogously, see below.) Writing 
as before k — \/k{k + 2) and inserting ([77]) , the Sk must satisfy (jH , which now 
takes the form 

Sk{t) + 3HSk{t) + {K^e~^"' + m^)Sk{t) = 0. (28) 
With the unit of time chosen so that H = 1, these equations simplify further as 
Skit) + iSkit) + (k^-^' + m^)Sk{t) = . (29) 



(IV) 

For each k = \/k{k + 2), the differential equation has two linearly inde- 
pendent solutions, given by 



y±^(0 = e-i*J±.(e-*«;) , (30) 



where the order of the Bessel functions is 



.^^^ (31) 

with the numerical value of the mass parameter m in the units employed, i.e. 
rn* = c^m/ (hH) when m is given in cgs units. The Bessel functions specialize 
in case that v = ±1/2 according to 

/ 2 / 2 

Ji{x)~\ — since J_i{x) — \ — cosx (32) 

2 V TTO; \ TTX 

for real x, and since this situation simplifies the numerical calculation of particle 
creation considerably, we will assume that m has been chosen such that v — 
±1/2. This corresponds to m* = v^, or m = ^/2hH/c^ « 3 x 10~^^kg in cgs 
units. This is surely a very small mass, but again, this value is chosen for the 
purpose of illustration so as to make the numerical calculations easier. With 
that choice of m and correspondingly = v^, and the identities p2p . one 
obtains for each fc e No a solution 

Su{t) = A{n)e'*- cos{e'* + B{K)e^* sin^^e"*/*) (33) 

to (|29p . where the coefficients 

, (2i + l)sinK cosK , 
A{k,) = ^ ^ + and (34) 

X — (2i + l)cosK sinK 
Bin) . -^-^ + ^ . (35) 

have been chosen such that the initial conditions 

Sk{0) = i, Sk{0) = ^, (36) 
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are fulfilled. This ensures that the Sk{t) defined in ([55]) fulfill the normalization 
condition (jM]) . and therefore give rise to a sequence {Sk} G 

(V) 

Now that we have chosen some {Sk} G •3'^, we will calculate the Bogoliubov co- 
efficients as described in (II) above. To this end, we must evaluate the integrals 

c^i\^) = l I dtMt){\Sk{t)\'+cul\Sk{t)\'- 



and 

which depend on the smearing functions fi{t) (i — 1,2) and Sk(t). After some 
algebra one can write the integrands as 



i^) = \ f dtm[Skitr+coiSk{tr) , 



\Skitr+ujt\Sk{t)\' = 



5 /3 15 \ 1 /3 15 

4 + T j + U ^ 2^=^ + U " " (37) 

K 5 \ 3 ^ 

4"4;^j^^ + 4;^^« 



and 

Sk{tf+LulSk{t)^ 



. 3 /3 9 3i \ f. A r 

^-4 + Tj^^+U~^^^y V^2j^^ (38) 
9 3i 



where the notation 



/ 3 i 
\Ak k 


K 

+ 4 






h 


— 


-2* cos [(2 - 


2e"*) k] 


h 


= e" 


3* sin [(2 - 


2e-*) k] 


h 


— 


-2* sin [(2- 


2e-*) k] 



h = e"3* cos [(2 - 2e-*) k] 

has been used. 
(VI) 

As weighting functions fi{t) we would have to use squares of C^-functions. 
However, since we are eventually evaluating the integrals c\ and numerically 
on a computer, it is justified to approximate (up to machine precision) weighting 
functions by Gaussians. Therefore, for our numerical calculations we choose as 
weighting functions fi{t) Gaussians localised at the times ti and characteristic 
smearing widths e and 6 as illustrated in Figure 2 (cf. also Figure 1): 

/i(t)^-e-(^) , /2(i) = Te-(^) . (39) 

e 
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A 





Figure 2. Illustration of the parameters ti,t2,e and 5 which characterize the 
weighting functions fi and f2- 

As said, this is a sufficiently good approximation for a weighting function as 
long as the smearing widths e and 6 are much smaller than the unit of time (i.e. 
much smaller than 1 with respect to our choice of units). The normalisation 
factors 1/e and 1/S have no impact on the jJ-'^\K) and A*^*'(k) since they will 
cancel. Nevertheless we wrote them down for one can better imagine the limit 
to the Dirac distribution when e and 5 tend towards zero. The time coordinate 
is not subject to any restrictions a priori. Therefore the integrals 



3i(e,ii) 



di/i(i)/i (j = l,...,6) 



have to be evaluated. In order to apply the identities 



' / — 2 
2 2 x/TT __EL. 

dxe~^ ^ sm\p{x + y)] = e sin[p2/] 

2 2 \/7r 

d.xe~^ ^ cos[p(a; + yf\ = e ^ cos\py] 



(40) 
(41) 



for the integrations, we carry out the Taylor expansion of the argument of the 
trigonometric functions to first order in t — ti. This approximation is justified 
if e <C 1, since then the integrand will only be important for |f — ti| <C 1. We 
obtain 

2k (1 - e-*) » 2k (1 - e-'i) + ^{t-ti) . (42) 



Multiplication of the two exponential functions yields 

e-(^)'e-«* = e-t^(*-(*i-t^'))%-«*i 



(43) 



where we used (ae)^ <C 1. We define a new integration variable r := t—{ti — ^e^) 
and rewrite the argument of the trigonometric function: 

2k ' 



2«(1 _ e-^*^) + _(t - f,) = _ _ _ - 1 + e*^) 



(44) 
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In our case we have 



p=^, ra«2K(l-e-*0- (45) 

Consequently, 



3i(e,ti) = ^ 




-4ti 














/ne 


-2ti 


















oxp(2ti) 


COS 


[2k(1- 


e' 




34(e,ti) = ^ 


\/Tre' 




^2^2 
oxp(2ti) 


COS 


[2k(1 - 


e 




35(e,ti) = ' 


\/ne' 




k2,2 
oxp(2ti) 


sin 


[2k(1 - 






J6(e,ii) = ^ 


\/ne 


-2tig- 


^2^2 
oxp(2ti) 


sin 


[2k(1- 


e~ 





(46) 



which completes our task of calculating the integrals Ci and for i — \. For 
i = 2 the same procedure as described above applies. We just have to perform 
the substitutions ti — > t2 and e ^ 5 and to carry over the corresponding as- 
sumption 5 <C 1. 

(VII) 

The explicit calculation of |/3fcp according to equation ist elementary but 
very cumbersome. Below, we contend ourselves with calculation of for 
several constellations of the parameters e, ii , (5 and t2 ■ The diagrams were gen- 
erated with the program "Mathematica" . The particle creation coefficients 
|/?(k)P = l/3/cp are plotted against n = \Jk{k + 2) continously, although ac- 
tually K is discrete in our case where E = 5'^. Furthermore, all times are 
represented in the natural unit tn = 1, 3787 x 10^°y. 

We briefly recall our assumptions made for our numerical calculations: 

• The scale factor is a{t) = e*. 

• The mass of the particle has numerical value — \/2 in the units used. 

• The test functions are approximated by Gaussians localised around ti,t2 
and have characteristic smearing widths e, S. 

• e,(5 < 1 

5.1 The case ti = t-2 

We consider the case ti = t2 and investigate the influence of e and 5 on the 
particle creation. The special case e = 5 is trivial: There is no particle creation 
since a;(l) = u}{2) 
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Figure 3. Particle creation rate for t\ = t2 = 0.1, 5 = 10 ^ for various values 
of e with e ~ 5 

As soon as e and 5 deviate minimally (e « S), small variations of e at constant 
S will affect the particle creation rate in the region k = e*i7r/(4e) ~ 

e*^7r/(4(5) (see Figure 3). Thus the magnitude of smearing widths determines 
which modes will be excited. As another case we consider a scale difference 
between e and S (e <C 5), cf. Figure 4. We remark that the maximum of the 
curve conforms to the larger smearing width (which is e in our case) 

according to the formula 

(47) 



4S 



Small variations of e influence |/3(k)P only around Hn 
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Figure 4. Particle creation rate for ti = t2 = 0.1, 5 = 10 ^ and various values 
of e for the case 6 
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Figure 5. Particle creation rate for e = 2^ = 2 x 10 ^ and logarithmically 
increasing values of ti = t2 

Now we look at the effects caused by increasing values of ti = ^2, keeping e 
and 6 fixed (Figure 5). With logarithmically growing preparation times, K^ax 
increases linearly, while |/3(tmax)P remains constant. Simultaneously the curve 
broadens. This is due to the properties of the solutions Skit) in an exponentially 
increasing Universe where the frequency of oscillations decays (x e~*. Thus 
the ratio of e and 5 (which remain constant) to the period of the oscillations 
decreases, leading to the same effect as shifting K,nax to the right which would 
be obtained by decreasing e and 5 at constant ^1,^2, cf. (jTZj) . Physically this 
corresponds to the redshift of a particle associated with the mode k, which loses 
energy. Thus at later times the same mode will be excited more easily. In the 
diagram this is represented by the fact that curves corresponding to later times 
arise from those at earlier times by a K-depending dilation to the right. 

Finally wc remark that the common feature of the case ti = t2 is |/3(0)p = 0. 

5.2 The case ti ^ 

We consider now the (more general) situation when ti and ^2 sue different. Due 
to the difference ^2 — ^1 there are phase differences in the trigonometric functions 
entering in |/3(k)P, which manifest themselves in oscillations of the particle 
creation curve as seen in Figure 6. Their frequencies grow with increasing 
time difference <2 — ^i- Due to these phase differences, |/3(0)p is no longer 
equal to zero. Thus, apart from the oscillation effects there is another new 
remarkable feature: the particle creation effect is biggest for small modes. When 
imposing much larger time differences ^2 — ^i, the particle creation curve grows 
and oscillates with higher frequency (Figure 7). 



20 



2.x 10"* 






1 1 1 








1 1 1 1 
- - - f2=0.1001 
-- r2=0.1004 


1.5 X 10"* 














r2 =0.1 007 ^ 


^ l.xlO"" 




" 1- . 












5.x 10"'' 




h. 


/' \- 


\ / \ 


























10000 20000 30000 40000 

K 



Figure 6. Particle creation rate for e = ^ = 10 ^ , ti = 0.1 and various values 
of t2 corresponding to a small time difference 

There is some similarity with findings in llj, where it is also argued that the 
particle number of a state prepared at early times can diminish when compared 
with a "vacuum state" at later time. However, the authors of TT use instanta- 
neous vacuum states as reference states, which is problematic since they aren't 
(locally) unitarily equivalent to Hadamard states [HI [H] • 
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Figure 7. Envelopes of the particle creation rate for e = S = 10~^, ti = 0.1 and 
various values of t2 corresponding to a big time difference. We have supressed 
the oscillations because they are of too high frequency to be properly represented. 

We have already remarked towards the end of Sec. 4 that in the case t2 > ^i, 
so that /i and /2 aren't proportional, it can't happen that |/3(k)P = for all 
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values of k. This is corroborated by the numerical results represented in Figures 
6 and 7. 

Another case of interest is the limit e — > and/or 5^0. In this case fi{t) 
and/or /2(t) would no longer be test functions, but delta distributions. Formally, 
the SLEs corresponding to the fi would, in these limits, tend to instantaneous 
vacuum states, which fail to have the Hadamard property. This implies that 
these limiting states are no longer unitarily equivalent to any other SLE. This 
in turn should be expressed in the divergence of the total particle number in the 
corresponding squeezed vacuum vector, which is equivalent to > 0{k,~^). 

In fact one can show (see 4^ ) that in our example we have for fixed i5 > in the 
limit e = 0, 

> ^ sin^ (2 (1 - + 0(«-3) . 

Comparing this with the criterion (PO)) . one can see that the limiting state for 
e — > isn't unitarily equivalent to the fixed SLE lo{2), and hence not a Hadamard 
state. 

6 Summary and outlook 

We have calculated the expected number |/3/cp of created particles per frequency 
mode k for an initial SLE at early cosmological time ti in a reference SLE at late 
cosmological time ^2 in a- closed, exponentially expanding Universe. A character- 
istic feature is that shows oscillatory behaviour with respect to variation 
of k. The envelope of decays in k more strongly than ^ , so that 

the oscillations are most significant for low frequency modes. The oscillatory 
behaviour increases with growing time-difference t2 — ti. 

A substantial drawback in discussing the possible physical significance of 
these findings is that they have been obtained under considerably simplifying 
assumptions which cannot really be considered as physically realistic. This refers 
mostly to having made the assumption of a spatially closed, exponentially ex- 
panding FRW spacetime and having fine-tuned the mass term m in the scalar 
field equation to an extremely small value so as to make the numerical calcula- 
tion more tractable. Above that, the quantized linear scalar field is a toy model 
to which there corresponds no observed particle. It is, then, much desirable to 
extend the investigation of this article to other quantum field modles such as 
the Dirac and electromagnetic fields, and for physically realistic cosmological 
scenarios and parameters. 

However, the basic features of our methods should carry over also to these 
more realistic models, with qualitatively analogous results. Therefore, the main 
purpose of this paper is to show that comological particle creation can be rigor- 
ously discussed, and detailed numerical results can be obtained, when employing 
states of low energy as reference states. This, in turn, also demonstrates the 
utility of the class of SLEs in the context of cosmological considerations. 
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